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ABSTRACT 

The Blandford-Znajek constraint equation for a stationary, axisymmetric black-hole force-free mag- 
netosphere is cast in a 3+1 absolute space and time formulation, following Komissarov (2004). We 
derive an analytic solution for fields and currents to the constraint equation in the far-field limit that 
satisfies the Znajek condition at the event horizon. This solution generalizes the Blandford-Znajek 
monopolc solution for a slowly rotating black hole to black holes with arbitrary angular momentum. 
Energy and angular momentum extraction through this solution occurs mostly along the equatorial 
plane. We also present a nonphysical, reverse jet-like solution. 

Subject headings: Black Hole Physics: Force Free Magnetospheres, Energy Extraction 



1. INTRODUCTION 

IPenrosd l|1969j) recognized the possibility to extract 
the spin energy of a black hole using particle de- 
cay in negative energy orbits within the ergosphere. 
Based on studies of force - free pulsar magnetospheres, 
iBlandford fc Znaiekl l|1977|) proposed that rotational en- 
ergy could be extracted through currents flowing in the 
black hole's magnetosphere. In this picture, strong elec- 
tric and magnetic f ields are induced by gravi to-MHD 
(GMHD) processes. IBlandford fc Znaiekl l)1977|) derived 
the equations for a stationary, axisymmetric force-free 
magnetosphere in curved spacetime, and reduced the 
set of equations to a central constraint equation relat- 
ing toroidal magnetic field H v to the charge density p 
and toroidal current density J v . They also found a per- 
turbative solution to the constraint equation valid in the 
limit a/M <C 1, where M is the black hole mass and a is 
t he angular momentum p er unit m ass. 

d Thorne fc MacDonaldl l|1982tt and 

acDonald fc Thornd tl 1982(1 developed this theory 
in a more intuitive "3 +1" formulation th at led to the 
membrane paradigm (T horne et al.l 119861) . where the 
equations of GMHD were written using the familiar 
electric and magnetic 3-vectors in absolute space whose 
time dependence is governed by Maxwell-type equations. 
IKomissarovl l(2004|) recently presented the essential equa- 
tions of this formulation in a form u seful for numerical 
studies, and helped resolve questions l|Punslv fc Coronitil 
119901 lPunslvl l 200 D relating to energy extraction in the 

membrane paradigm. 

I ■ f equations presented by iKomissaro 2 il200l pro- 
vide a useful starting point to search for analytic so- 
lutions. Here we us e these equations to reder ive the 
constraint equation of lBlandford fc Znaiekl ltl977l) in the 
3+1 form. This brings forth a clear understanding of 
the nature of the poloidal functions defining the currents 
and fields. We have discovered an analytical solution 
valid for arbitrary angul ar momentum that r e duces to 
the monopole solution of Bland ford fc Znaiekl (jl977fl in 

1 Troy University, Troy, Alabama, 36082 

2 E. O. Hulburt Center for Space Research, Code 7653, Naval 
Research Laboratory, Washington, DC 20375-5352 

3 On sabbatical leave from Troy University. 



the lim i t of a/ M <C 1. This solution, which satisfies the 
iZnaiekl (|1977f) regularity condition, permits energy ex- 
traction preferentially along the equatorial direction of 
the Kerr black hole. 

The modified Maxwell's equations in curved space-time 
are given in Section 2, and the equations for a force-free 
magnetosphere are presented in Section 3. In Section 4, 
we construct the form of fields and currents for a given 
poloidal function Q. The governing constraint equation 
for this function is given in Section 5. Solutions to this 
equation are derived in Section 6, and we summarize in 
Section 7. 

2. ELECTRODYNAMICS IN ABSOLUTE SPACE 

While Maxwell's equations preserve all of its elegance 
in a covariant formalism on a four dimensional mani- 
fold, it distracts from some of the simple (far-field) so- 
lutions that it might permit. With this is mind, we 
briefly state the essential equations of electrodynamics 
in an absolute thre e dime nsional space. The recent pa- 
per by IKomissarovl l|2004|) explains how these equations 
are derived. 

The construction of absolute space is facilitated by not- 
ing that an arbitrary spacetime metric can be written in 
the form 



ds 2 = (f3 2 - a 2 )dt 2 + 2f3 i dx i dt + ^ ij dx i dx j . (1) 

The functions x % serve as coordinates for our spacelike 
hypersurfaces defined by constant values of t. Consider 
one such hyperspace £ defined by the region t = 0. We 
can think of electric and magnetic fields (E and B) as 
objects existing in our absolute space E. The time evo- 
lution equations for E and B in the presence of a charge 
density p and electric current density vector J in our 
absolute (curved) space endowed with a metric 7^ are 
given by the following set of Maxwell's equations: 



V ■ B = 0, 

d t B + V x E = 0, 
and their inhomogencous counterparts, 



(2) 



(3) 
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V • D = p, (4) 



-d t D + VxH = J. (5) 

It is important to remember that E, B, D, H, and J are 
vectors in our three dimensional absolute space £, and 
in general are time dependent. Also, V is the covariant 
derivative induced by the metric jij on S. As usual, the 
curl of a vector field is defined by the expression 



(V x Af = e'^VjAk, 



(6) 



where e l ^ k is the completely antisymmetric pseudotcnsor 



such that e 123 



, and 7 = Det(7y). It is easily seen 



that Maxwell's equations imply the continuity equation 



dtp + V • J = 0. 



(7) 



Unlike its flat space counterparts, even in regions of neg- 
ligible electric and magnetic susceptibilities, E ^ D, and 
B 7^ H . Indeed, it can be shown that they instead satisfy 
the consitutive relations 



and 



E = aD + f3 x B, 
H = aB - f3 x D. 



(8) 
(9) 



Of interest are spacetimes admitting Killing fields cor- 
responding to axial symmetry (m) and stationarity. Con- 
sequently, Noether's theorem imply energy and angular 
momentum conservation laws. They can be stated in the 
form 



and 



Here, 



d t e + V- S= -(E- J), 



d t l + V ■ L = -{pE + J x B)-m. 



e= ~(E ■ D + B ■ H) 



is the volume density of energy, and 
I = (D x B) ■ ni 
is the density of angular momentum, 

S = E x H 
is the flux of energy, and 

1 



(10) 
(11) 
(12) 

(13) 
(14) 



In our absolute space framework, the force free condi- 
tion reduces to 



and 



(16) 
(17) 

These restrictions, along with Maxwell's equations and 
Eqs. © and ©, imply that 



E-J = Q, 



pE + J x B = 0. 



and 



E T = 0, 
E P - B P = 0. 



(18) 
(19) 



The poloidal and toroidal components (Ap, and At) of 
a vector field are defined such that A = Ap + At, where 
A P = A r d r + A e d e and A T = A^d v . Eqs. (JTSJl and ((THJ 



imply that there exists a vector u> 



Qd v such that 



E 



- U! X B, 



(20) 



From the vanishing of the curl of E under the stationarity 
condition (Eq. ©), one finds that 



B ■ V£! = 0. 
It can also be shown that 

B ■ VH V = 0. 



(21) 
(22) 



4. EXPLICIT EXPRESSIONS FOR FIELDS AND CURRENTS 

To simplify calculations, we shall assume that the spa- 
tial coordinates are orthogonal, and that the shift vector 
is purely toroidal, i.e., (5 = (0,0, (3^), The Kerr solu- 
tion written in Boyer-Lindquist (though not Kerr-Schild) 
coordinates can be written in this form. 

Surfaces of constant O are referred to as poloidal sur- 
faces (not to be confused with poloidal components of a 
vector). From Eq. (|21|l it is clear that B is tangent to 
poloidal surfaces. Since ft does not have any ip depen- 
dence, and since Eq. Q21JI has nothing to say about the 
toroidal component of £>, it is clear that Bp will entertain 
solutions of the type 



Bp — —=(—Qed r 

V7 



VL^ r dg) 



(23) 



where, for the moment, A is an arbitrary function. This 
must be so because in the two dimensional subspace given 
by Q — const, there is a unique vector (modulo magni- 
tude) that is perpendicular to Vf2. The condition that 
B is divergence free means that A satisfies 



L = -(E ■ m)D - (H ■ m)B + -{E ■ D + B ■ H)m (15) 
is the flux of angular momentum. 

3. STATIONARY, AXISYMMETRIC FORCE FREE 
MAGNETOSPHERES 

The condition that the magnetosphere is force free 
brings about enough structure into Maxwell's equations 
to enable the introduction of a streaming function that 
will help us visualize the field structure in geometric 
terms. It is traditional to use spheroidal spatial coor- 
dinates given by x l = (r,6,<p) such that m = d v . As- From ©, we see that 
sumptions of stationarity and axissymmetry imply that 
d v g^u = = dtg^v 



A r ft a = A ail. 



(24) 



Consequently, A is a poloidal function (a function that is 
constant on poloi dal surfaces). In the notati on of the 
original paper by IBlandford fc Znaiekl lll977l) . AdQ = 
—dA v . The electric field is immediately calculated from 
Eq. ()20J) and, as expected, is the gradient of a scalar 
function: 



E P = Ad{fl 2 /2) 



A 



Andn. 



d = d p = -(n + f3 v )dn. 

a 



(25) 



(26) 
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Similarly, the expression for Hp can be calculated from 
Eq. ©, giving 



Hp = (a 2 - 1 - — 
a 



(27) 



The electric charge is determined by the divergence of 
the Dp, (Eq. @). Explicitly, 



pE ■ E p + (J x B) ■ E P = 
P E 2 + {{J P + J T ) x (Bp + B T )) • E P = 
P E 2 + ((Jp x B T ) + [Jt x Bp)) -E P = 0, (33) 

since Jp is parallel to Bp. With the help of the following 
relations. 



Jp = 



(28) 



}_dHf 
77 dVL 



(Slfidr — $l, r de 



dH, 



B f 



The toroidal component of the electric current density 
vector can be obtained from the derivatives of compo- 
nents of Hp : 



A 



E r Bn — Eh B r 



nA 2 

7^ 



(lee{n r ) 2 + lrr (n e ) 2 ), 



and 



r! 2 A 2 



E 



■(j ee (n, r ) 2 + lrr (tt fi ) 2 ), 



^J v - H g . r - H r .e - d r [——(a - 13 - l3 lp n)j eg ri ir }+ Eq ^ re d UC es to the manageable form 



do 



A 



{a 2 -p 2 - (3 v n) lrr n, 



(29) 



It clear from the above discussion that the poloidal fields 
and, consequently, the toroidal current J v are uniquely 
described by the poloidal functions and A. On the 
other hand, the toroidal fields and the poloidal currents 
can be determined from the poloidal function H v . In 
particular, from Eq. 0, it is clear that H v = aB v . 
Maxwell's equation (Eq. JSJ) implies that 



y/jJp = H^fid r — H vr de. 



(30) 

Thus we see that fields and currents separate into two 
distinct categories: objects that are determined by Vl 
and A, and those that are determined by H v . Outside of 
the fact that H v is a poloidal function (by definition, Q 
is), it is not yet clear as to how these two functions are 
dynamically related. This issue will be cleared up in the 
following section. 

5. THE CONSTRAINT EQUATION 

The expressions for the fields and currents given in 
the previous section naturally satisfies Eq. (|16|l . Since 
the toroidal component of the electric field vanishes, it 
is easily checked that, from Eq. Ijl7jl . (J x B) v = (as 
shown below in Eq. H32|0 . Thus the only remaining re- 
quirements for a force-free solution is 



pE P + (J x B) P = 0. 



(31) 



The implication of the above equation is most easily un- 
derstood by projecting the equation onto Ep, B p , which 
serve as a basis vectors for poloidal vector fields. The 
above equation yields no constraint when projected onto 
Bp, i.e., 

pE ■ B p + (J x B) ■ Bp = (J x B) ■ (B - B T ) = 

- B*(J x B) = -5"Ai(I ft( ll |r - Hp rQ e) = 0. 

(32) 

Projecting Eq. I|17|) onto Ep gives 



1 dH 2 

2A dn 



a(pQ-f, 



'J U2 



(34) 



(35) 



This is the final and only constraint equation. If Q and 
A are picked such that the right hand side of the above 
equation is a poloidal function, then H v continues to be 
poloidal function. The poloidal functions J7, A, and H v 
then uniquely determines all currents and fields. It is 
important to realize that f2 is not to be thought of as a 
potential: physically relevant quantities like the electric 
field depend on fl directly, and are not invariant transfor- 
mations of the type O — s- f7 + const. The charge density 
p and the toroidal current J v are functions of and A 
(see Eqs. (EHJ) and (PI)- 

6. ASYMPTOTIC SOLUTIONS AND ENERGY 
EXTRACTION 

Inserting Eqs. (JSHJ and (J2HJ into Eq. our con- 

straint equation gives 



1 dH 2 
2A dfl 



22s£[n*(-A 



<V7 



;(7^fi + /3 v ,)j ee n :r )+ 



A 



0r( 



A 



A 



(S3 2 -a 2 



(P 2 



a 



■ f3 v Sl)^ 6 ett,r)- 



(36) 



Therefore, Eq. 



(HSJ is equivalent to Eq. (3.14) of 
IBlandford fc Znaiekl l|1977|) written in the 3+1 formal- 
ism. 

While searching for solutions for the fields and currents 
that might permit extraction of energy and angular mo- 
mentum from a rotating black hole, it is advantageous to 
observe that 



d 2 E 
dAdt 

d 2 C 
dAdt 



= S r 



U 



Hr 



— — H^p£lB r 



-H ip B r 



Hr 



(37) 
(38) 
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as can be see from Eqs. (|14fl and l|15(l . Here £ and C are 
the total energy and angular momentum, respectively, 
extracted from the black hole. 

For definiteness, we shall consider the magnetosphere 
of a Kerr black hole in Boyer-Lindquist coordinates. For 
finite rates of energy and angular momentum extraction, 
it is clear from the above two equations that for r ^> 
M, CI — > Cl(9). With this in mind, we seek solutions 
to the constraint equation of the type CI = £1(0) for all 
values of r. This means that all poloidal functions are 
functions of alone, since all poloidal functions are of 
"zeroeth" order in r. Due to the inherent complexity 
of the constraint equation, we shall further consider Eq. 
(|36l) in the far field limit. To order (1/r 3 ) for strictly 
^-dependent functions, Eq. I|36|l takes the form: 



1 



dH 2 



2/(0) d9 



- ft sin — (/ft sin ( 
dv 



-^[-o^sm^C/fism*) + ^(^)] 
H dti dO sinO 



+2M— [aft-^-(/sin0(l 
r J dt) 



an sin 2 0))- 



d0 V sin0 V 



aft sin 2 0))], 



(39) 



where M and a are the mass and the angular momen- 
tum per unit mass of the black hole, respectively, and 
f{9) = — Mlfi = A Vy g. For a consistent formulation of 
the theory of axisymmetric, stationary, force-free mag- 
netospheres, the above equation implies that if H v is to 
be a poloidal function of alone, then the terms propor- 
tional to the inverse powers of r must vanish identically 
for choices of / and ft. 

General solutions to order 1/r 2 can be considered by 
ignoring the 1 /r 3 term in the right hand side of Eq. (I39() . 
Here we require that ft and / satisfy the relation 



a 2 ftsin 2 0— (/ftsin0) = — (- 



dO 



dQ sin ( 



(40) 



To solve this, let g = /ft sin0 and h = (ft sin 2 0)" 1 . With 
these definitions, Eq. I|4U|I becomes 



a 2 d d . 



(41) 



Integrating the above equation results in the relation 

Ciftsin 2 



fJ 



V\ a 2 ~h 2 \ 
Consequently, for an arbitrarily chosen ft, 



/ 



Ci sin t 



(aft sin 2 0) 2 - 1 



(42) 



(43) 



would make the 1/r 2 term in Eq. I|39[) vanish. We can 
therefore successfully obtain the following function for 
H v . Explicitly, 



The choice of ft is determined by the Znajek regular- 
ity condition applied at the event horizon (r+ = M + 
V M 2 — a 2 ) so as to make B v finite in the well-behaved 
(even ne ar the event horizon) Kerr-S child coordinate sys- 
tem (seeEHiiiS Jll23);|KHmisiar23 (EH)- The Znajek 
condition can be written as 



a. 



sin 



-(2r+Mft- a)B r 



i0 



(2r+Mft - a)/, 



( 45 ) 

where the subscript + indicates that the relevant quan- 
tities are to be evaluated at the event horizon and p 2 ^ = 
r\ + a 2 cos 2 0. From Eqs. and (|43|l . we see that 



\(4:r 2 + M 2 - p%)Cl 2 -Ar + MaCl+a 2 ]f . (46) 



We shall consider the solution for H v such that Hq = 
(it is easily seen that when H ^ 0, the resulting solu- 
tion does not permit a finite rate of energy extraction, 
and this type of situation will be dealt with in Subsec- 
tion 6.2). This is possible if and only if the quantity in 
the square brackets in the above equation vanishes iden- 
tically. Solving the resulting quadratic equation for ft, 
we find two solutions, namely 



Q 4 



2Mr+ + p\ 



and 



ft_ 



2Mr, 



p\ a sin 2 I 



(47) 



(48) 



From Ea. l|43|l and the definition of /, we see that the 
only non-vanishing poloidal component of the magnetic 
field is given by 



B' 4 



1 1 B sm6{2Mr + ± p 2 + ) 

V7 (a sin 0) 4 - (2Mr+ ± p\ f \ 



(49) 

where we have relabeled C\ as B$. It is clear that ft_ is 
an unphysical solution since B r as given above is unde- 
fined everywhere. 

6.1. The ft + Solution 
In this case, the non- vanishing components of the fields 



B r = 



1 B sin \/oTTh 



Eg = -^ci + B r 

H v = -^Cl+B r sm6, 



(50) 



where ft# — a/2Mr + is the angular velocity of the event 
horizon. When a <C M 

B r , — > Bo sin , and 



dH 2 d 
-^2/ftsin0-(/ftsin 



de 

(/ft sin 0) 2 



(/ft sin t 



8M 2 



(51) 



This is precisely the iBlandford & Znaiekl l)1977f) 
(44) monopole solution ijKomissaro vl 120041) . Therefore, the 
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solutions for the fields and currents corresponding to 
£1 = f2+ generalizes the Blandford-Znajek monopolc so- 
lution to accommodate the case of a black hole for all 
values of a 2 < M 2 . 

A parallel approach to the study of the force-free mag- 
netosphere has been developed v ia the Grad-Shafranov 
equati on (see, e.g., Eq. (6.4) of iMacDonald fc Thornd 
(1982)). In our notation, the G rad-Shafranov equation 
takes the form (Uzdcnskv 2005) 



y [ oV^ 1 (Sl + + ^) 2 lw )] 



7vv 



dij: 



W) 5 



I dl 

aj vv dip 



0. 



(52) 



Here, I — H v and tp — A v . By straightforward substi- 
tution and evaluation of the various terms in eq. (|52l) . 
it is not difficult to see that our solution satisfies the 
Grad-Shafranov equation to order 1/r 2 . 

From Eqs. I|50|) and (|37fl . the angular dependence of 
energy extraction can be calculated. In the limit r ^S> M, 
the result is 



d 2 £ 
dAdt 



clQh ( B$ 2 sin 



(53) 



From the above equation, it is clear that most of the 
energy extraction happens along the equatorial plane. 
The total rate of energy extraction can be obtained by 
integrating the above result, giving 

= _o arctan— - — . 54 
at ar+ r+ 2M 

In similar fashion, we see by integrating Eq. I|38l) that 



dC 
~dt 



2n 



Bin 



1 d£ 
Qh dt 



(55) 



As a result of energy and angular momentum extrac- 
tion from the black hole, the mass and the total angular 
momentum (J = aM) of the black hole changes by the 
amount 

SM d£ SJ dC 

and = — — , (56) 



St 



dt 



St 



dt ' 



respectively. From Eq. (|55|l and the above definitions, it 
clear that 

SJ 2tt 2o 1 SM 
-^ + Y B 2 n H = — — . (57) 



n H st 



Ther efore we get the familiar inequality ijChristodouloul 
Il970l) 

SJ 1 SM 

H * t£-6t> (58) 

which ensures that the irreducible mass of the black hole 
is non-decreasing if the black hole evolves along a Kerr 
sequence in a reversible way. This process therefore can- 
not lead to the formation of a naked singularity. 

6.2. A Jet-Type Solution 

It is easily seen that SI = S7_ removes all the in- 
dependence in the right-hand side of Eq. (|39"|) to order 
1/r 3 . As shown by Eqs. igBJ and (15jl. SI = Sl_ is not 
a physical solution for the condition Hq = 0. We now 
let Hq 7^ 0, and impose the Znajek condition, Eq. (|45l) . 



for this case. Because our solutions involve both Q + and 
Q_, the results continue to be valid only to order r~ 2 . 
From Ea. (|46|l we see that 



±H 2 P % n+n 2 _ 



(59) 



a (fi-n + )(n-n_)' 

Here the ± factor is to ensure that f 2 > 0. Similarly we 
find from Eq. gHJ) that 



(60) 



a 2 sm 2 e 1 (n-n_ )(fi + n_) |' 

Equating the right-hand sides of the last two equations, 
we see that 



±H 2 p\ 



sin 9(2Mr + + p z + ) 
It is important to remember that any SI satisfying the 
above equation is consistent with Eq. (|39|1 (to order 1/r 2 ) 
and with Eq. (|45p. The above equation has the unique 
solution 

Au + fln- (62) 



A — B 



where 



A_ 



-1, 
-1, 



if n p - si+ > 

otherwise 



and 



B 



if n p + n_ > 

otherwise. 



(63) 



All other poloidal fields quantities are now uniquely de- 
termined by noting that / is given by Eq. (|60|l . It is 
important to see if we can indeed satisfy the above con- 
ditions. A quick calculation shows that 



Sip — Sl_j_ — 



o_ = 



A — B 

i(si+ 4 



and 



(64) 



-H%p\n+n- is 



B 2 a 2 sin 4 0] 



A-B 

Therefore the choice A — —B 2 and B 
a valid one. We shall pick this choice for the remainder 
of the paper. Consequently, we have 

si p = -n + n_ } H ° p + — -»~ — ; J . i(». r ,i 

P [H 2 p\n+ + B 2 a sin 2 9} V ' 

Note that as 9 — > and tt, H p — » — Sl_ — > —00 . The 
form of / is determined by Eq. I|60[) . and upon substitu- 
tion of the explicit form of SI (= Q p ), we find that in the 

limit as 9 — > and tt, f — > ±i?o ya/ftij/2. 

The expression for the rate of total energy extraction 
is given by 

d2£ - o 1 f ~ 1 f 2 o ( 2Mr +^ - fl ) f fif n 

pj • (66) 

As 9 — > and tt, 
d 2 £ 



■1 ff^a 

^Isi^" 



si_(si_ + si ff ) 



(67) 



A solution of this type has the following features: Energy 
extraction is less than zero near the poles, i.e., energy is 
being fed into the system, indicating a reverse jet type 
situation. Also, the total rate of energy and angular 
momentum "insertion" is not calculable since the above 
integral is divergent along the poles. This solution is 
therefore unphysical. 
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7. DISCUSSION 

Bas ed on the 3+1 equations as written bv lKomissarovl 
( 2004) , we have rederived the constraint equation relat- 
ing the toroidal magnetic field to the charge and current 
densities in a force- free magnetosphere around a spinning 
black hole. Known solutions to the constraint equation 
for the force-free magnetosphere include the monopole 
and th e parabolic solutions obtain ed in the orginal pa- 
per bvlBlandford fc Znaiekl l)1977j) . and the solution by 
iBeskin et alJ l|1992|) for a black hole surrounded by a 
magnetized, conducting accretion disk. We have dis- 
covered a solution to the constraint equation that gen- 
eralizes the "monopole sol ution" originally derived by 
Blandford fc Znaiekl (JT977). This solution satisifes the 
Znai eld 111 9771) regularity condition at the event horizon, 
even in the limit o/M<l (contrary to the statement of 
B landford and Znaje k) . 

Komissarov (2001) has used a time-dependent numeri- 
cal simulation to calculate the electromagnetic extraction 
of energy for a monopole magnetic field at different val- 
ues of a/M. Our value of SI + /Qh at 6 = 0.5 ranges from 
0.5 to 0.58 when a/M varies from 0.1 to 0.9, in compar- 
ison with the numerical value of 0.52 for a/M = 0.9 at 
r = 10. The value of H v for our £1 + solution is « 25% 
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using the techniques developed in this paper. 
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